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We report the creation of an ultracold neutral plasma by 
photoionization of laser-cooled xenon atoms. The charge car- 
rier density is as high as 2 x 10® cm~"^, and the temperatures 
of electrons and ions are as low as 100 mK and 10 ^K, respec- 
tively. Plasma behavior is evident in the trapping of electrons 
by the positive ion cloud when the Debye screening length 
becomes smaller than the size of the sample. We produce 
plasmas with parameters such that both electrons and ions 
are strongly coupled. 



The study of ionized gases in neutral plasma physics 
spans temperatures ranging from 10^^ K in the magneto- 
sphere of a pulsar to 300 K in the earth's ionosphere [Q . 
At lower temperatures the properties of plasmas are ex- 
pected to differ significantly. For instance, three-body 
recombination which is prevalent in high temperature 
plasmas, should be suppressed If the thermal energy 
of the particles is less than the Coulomb interaction en- 
ergy, the plasma becomes strongly coupled, and the usual 
hydrodynamic equations of motion and collective mode 
dispersion relations are no longer valid |^ . Strongly cou- 
pled plasmas are difficult to produce in the laboratory 
and only a handful of examples exist Q, but such plas- 
mas do occur naturally in astrophysical systems. 

In this work we create an ultracold neutral plasma with 
an electron temperature as low as Tg = 100 mK, an ion 
temperature as low as Ti — 10 /iK, and densities as high 
as 71 = 2x10^ cm"'^. We obtain this novel plasma by pho- 
toionization of laser-cooled xenon atoms. Within the ex- 
perimentally accessible ranges of temperatures and den- 
sities both components can be simultaneously strongly 
coupled. A simple model describes the evolution of the 
plasma in terms of the competition between the kinetic 
energy of the electrons and the Coulomb attraction be- 
tween electrons and ions. A numerical calculation accu- 
rately reproduces the data. 

Photoionization and laser-cooling have been used be- 
fore in plasma experiments. Photoionization in a 600 K 
Cs vapor cell produced a plasma with Te > 2000 K [^, 
and a strongly coupled non-neutral plasma was created 
by laser-cooling magnetically trapped Be+ ions 

A plasma is often defined as an ionized gas in which 
the charged particles exhibit collective effects [Q. The 
length scale which divides individual particle behavior 
and collective behavior is the Debye screening length Ad . 
It is the distance over which an electric field is screened by 



redistribution of electrons in the plasma, and is given by 
\d = \J ioksTe/e'^n. Here, eg is the electric permittivity 
of vacuum, fcs is the Boltzmann constant, and e is the 
elementary charge. An ionized gas is not a plasma unless 
the Debye length is smaller than the size of the system 
[Q. In our experiment, the Debye length can be as low 
as 500 nm, while the size of the sample is ct « 200 /im. 
The condition A^i < a for creating a plasma is thus easily 
fulfilled. 

The atomic system we use is metastable xenon in the 
6s [3/2)2 state. This state has a lifetime of 43 s 1^ and 
can be treated as the ground state for laser-cooling on 
the transition at 882 nm to the 6p[5/2]3 state [Q. The 
metastable atoms are produced in a discharge and sub- 
sequently decelerated using the Zeeman slowing tech- 
nique. The atoms are then collected in a magneto-optical 
trap and further cooled with optical molasses to approx- 
imately 10 /iK. We characterize the cold neutral atoms 
by optical absorption imaging [|l0|. This measurement 
provides the density and size of the atomic cloud and the 
number of atoms in the sample. Typically we prepare a 
few million atoms at a density of « 2 x 10^" atoms/cm^. 
Their spatial distribution is Gaussian with a rms radius 
a « 200 fim. 

We partially ionize the cold atom sample via two pho- 
ton excitation. A pulse of light from the cooling laser 
at 882 nm populates the 6p[5/2]3 level. Green pho- 
tons (A = 514 nm) from a pulsed dye laser, pumped 
by a frequency-tripled pulsed Nd:YAG laser, then excite 
atoms to states at or above the ionization potential. 

The energy difference, Ai?, between the photon energy 
and the ionization potential is distributed between the 
electrons and ions. Because of the large ion to electron 
mass ratio, all except 4 x 10~^Aii^ is given to the elec- 
trons. Equipartition of energy between ions and electrons 
requires tens of ms jl^]. We vary AE/ks in a controlled 
manner between 0.1 — 1000 K by changing the green laser 
frequency. The bandwidth of the laser of 0.07 cm^^ sets 
the lower limit. 

By adjusting the pulse energy of the green laser, we 
control the number of atoms photoionized. For the high- 
est energy available, 1 mJ in a 10 ns pulse, we can pro- 
duce up to 2 x 10^ ions, which corresponds to a peak 
density of 2 x 10^ cm^^. The number of atoms photoion- 
ized varies linearly with the laser intensity. Although the 
ionization fraction is low (< 10%), the charged particles 
show no evidence of interactions with the neutral atoms. 
This is to be expected because the mean free path for 
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neutral-charged particle collisions is much greater than 
the sample size ||l2|-|l|l. 

For detection of charged particles, an external elec- 
tric field is applied to direct ions towards a microchan- 
nel plate detector and electrons towards a single channel 
electron multiplier. The efficiencies are 50% for ions p5| ] 
and 85% for electrons. The magnitude of the apphed 
electric field is calibrated through field ionization of Ry- 
dberg atoms p6[ |. 




Time [fis] 

FIG. 1. Electron signals recorded for four different pulse 
energies of the green laser, i.e. different densities of charged 
particles (10^ — 10^ cm~^). The uppermost curve corresponds 
to the lowest energy. The photoionization occurs sd, t — 0. 
The initial kinetic energy of the electrons is AE/ks = 0.6 K. 
The data shown is an average over 20 cycles of the experiment. 
Also shown is the magnitude of the applied electric field. 

In each cycle of the experiment the atoms are first laser 
cooled and an electric field of approximately 0.005 V/cm 
is applied. The atoms are then photoionized, and after 
about 500 ns of time of flight a pulse of electrons arrives at 
the detector (see Fig. |^). If the green laser energy is high 
enough, the first peak develops a tail, and a second peak 
appears when the electric field is linearly increased a few 
microseconds later. On this time scale the ions are es- 
sentially stationary. About 300 /iS after the electric field 
ramp is applied, they are detected on the microchannel 
plates. 

A simple model (Fig. |2|) explains the experimental 
data. The charge distribution is everywhere neutral im- 
mediately after photoionization. Due to the initial kinetic 
energy of the electrons (« AE), the electron cloud begins 
expanding. The resulting local charge imbalance creates 
an internal electric field which produces a Coulomb po- 
tential energy well for electrons. If the well never be- 
comes deeper than the initial kinetic energy, all the elec- 
trons escape. This corresponds to the uppermost curve 
(i.e. lowest laser intensity) in Fig. ^. If enough atoms 
are photoionized, however, only an outer shell of elec- 
trons escapes, and the well becomes deep enough to trap 
the rest. Electrons in the well redistribute their energy 
through collisions within 10 — 100 ns As charges are 



promoted to energies above the trap depth they leave 
the well. This explains the tail of the first peak in the 
electron signal. During this process of evaporation the 
potential well depth increases. Evaporation eventually 
slows and remaining electrons are held until an applied 
electric field overcomes the trapping potential. They ap- 
pear as the second peak in Fig. |l]. 

This description suggests that for a given AE there 
is a threshold number of positive ions required for trap- 
ping electrons. The data show such behavior in a plot 
of the fraction of electrons trapped versus the number of 
photoions produced (Fig. ||a). As AE increases, more 
positive charges are required to produce the trapping ef- 
fect. 

At the trapping threshold, after all the electrons have 
left, the potential well depth equals the initial kinetic 
energy of the electrons. From this relation one can 
calculate the number of positive ions at the threshold, 
N* = AE/Uq. Here, Uq = ^2/^ eV47reoCT, and a is the 
rms radius of the Gaussian spatial distribution of positive 
ions. 
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FIG. 2. Schematic of the potential energy seen by a test 
electron when enough atoms are photoionized to result in 
trapping of electrons. Photoionization occurs at to = and 
the sample is everywhere neutral. Because of the kinetic en- 
ergy imparted by the laser, some electrons leave and a charge 
imbalance develops. At ti « 10 ns the resulting potential 
well equals the initial kinetic energy, trapping the remaining 
electrons. Due to Debye screening the bottom of the well is 
fiat. As electrons in the well thermalize, evaporation occurs. 
The well depth increases and the electrons cool slightly. By 
t2 ~ I fis evaporation essentially stops. The dashed line indi- 
cates the average kinetic energy of the electrons. 

The data agree with this simple calculation. For a wide 
range of electron kinetic energies, the onset of trapping 
occurs at = N* , as shown in Fig. ||b. Scaling the 
number of ions produced by N* shows that all data fall 
on a universal curve p^ . A numerical simulation which 
approximates the initial velocity, v = \J2AE jm^ as di- 
rected radially outward and integrates the equations of 
motion for the electrons, reproduces this behavior. 
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The model described above (Fig. ^ implies that the 
temperature of the electrons is Te^AE/kB- This is con- 
firmed by the numerical simulation. In principle, the 
energy distribution of the trapped electrons can also be 
determined from the shape of the second peak in Fig. 
|l|. However, the analysis is complicated because the trap 
depth increases as electrons are removed. Also, rether- 
malization is fast and the temperature changes on the 
timescale of the electric field ramp. 

The initial ion temperature is easily estimated. For 
excitation close to the ionization threshold, the energy 
imparted to the ions from photoionization is negligible 
compared to the kinetic energy of the atoms. Therefore 
the minimum initial temperature is 10 /xK. For large AE 
the temperature approaches 4 x 10~^AE/kB, which is 
4mK for AE — 1000 K. Although the equilibration time 
is on the order of tens of ms, collisions with the energetic 
electrons are expected to approximately double the ion 
temperature within a fis. 

The threshold condition N = N* is mathematically 
equivalent to Xd = a. In this context, one can interpret 
Ad as the displacement of electrons from their equilib- 
rium positions when their energy in the local internal 
electric field in the plasma equals their kinetic energy 
0. If A/) > (7, the electrons are free to escape to infinity. 
li Xd < electrons are trapped in the ion cloud by the 
internal field and a plasma is formed. 

After the untrapped fraction of electrons has escaped, 
the cloud as a whole is no longer strictly neutral. But as 
mentioned above, electrons escape most easily from the 
edges of the spatial distribution, and for N > N* the 
center of the cloud is well described as a neutral plasma. 
This behavior is also seen in the numerical simulation. 

The only significant effect of the residual charge im- 



balance is a Coulomb expansion of the cloud that occurs 
on a long time scale of many microseconds. This limits 
the time available for studying the highest density con- 
ditions. The expansion also decreases the potential well 
depth, allowing formerly trapped electrons to escape. In 
a plasma with 5000 ions and 10% charge imbalance, half 
of the initially trapped electrons escape in about 100 fis. 
The expansion is slowed compared to what would be ob- 
served for a bare cloud of positive charges, however. A 
bare cloud of 5000 ions, initially with cr — 200 /xm, ex- 
pands to twice its radius in a few microseconds, reducing 
the well depth by a factor of two. 

Phenomena similar to the electron trapping observed 
here are seen in traditional plasmas. For instance, recom- 
bination can often occur at containment walls and leads 
to net charge diffusion from the center of the plasma. 
The mobility of the electrons is larger than that of the 
ions, but their motion is retarded by local internal elec- 
tric fields which develop from any charge imbalance. This 
leads to ambipolar diffusion [Q in which electrons and 
ions migrate at equal rates. 

In our ultracold plasma the thermal energy of the 
charged particles can be less than the Coulomb in- 
teraction energy between nearest neighbors, making it 
strongly coupled. The situation is characterized quan- 
titatively by the electron and ion Coulomb coupling 
parameters Q Fg = (e^/47reo a)/fcs?e and F^ = 
e-^/^^FeTe/T, where a = (47rn/3)-i/3 is the Wigner- 
Seitz radius. The exponential term in the expression for 
Ti is due to the shielding of the ion-ion interaction by 
electrons [Q. When F > 1, many-body spatial correla- 
tions J2G] exist and phase transitions such as crystalliza- 
tion ]21| ] of the sample may occur. Systems with Fg > 1 
are sometimes called non-Debye plasmas because Ad < a. 
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FIG. 3. (a) The fraction of electrons trapped is plotted versus the number of photoions created. Each curve corresponds to 
a different green laser frequency. The corresponding initial energies of the electrons are displayed in the legend. As the energy 
increases, more positive ions are required to trap electrons, (b) Same as (a) but the number of photoions is scaled by A'^*. The 
threshold for trapping is given by A'^ = A'^*. The line is the result of a numerical simulation. There is a scale uncertainty of 
about 10% in determining the fraction of electrons trapped. 
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For the densities and temperatures accessed in the ex- 
periment, we can prepare a plasma in which both elec- 
trons and ions are initially strongly coupled: Fe = 10 and 
Ti = 1000. To our knowledge such a system has never 
been created before. 

This novel plasma is well suited for a wide range of 
experiments. Plasma oscillations [0,^, which have a 
frequency fp — \/r«?/mio/27r of up to 400 MHz, can 
be used to probe the density distribution of the system. 
Magnetic confinement may greatly extend the plasma 
lifetime, and because of the low sample temperature, the 
required field should be small. The thermalization and 
evaporative cooling of electrons, and the temperature of 
the ions require further study. 

We can also look for three-body electron-ion recombi- 
nation. At higher temperatures, the rate for this process 
scales as T^^/^ and an extrapolation to the present 
experimental conditions yields a recombination time of 
nanoseconds (for Te — IK and n = 2 x lO^cm"^). The 
long lifetime we observe 100 /is) is the first clear indi- 
cation that this theory, as well as an extension to T « 1 K 
breaks down for the temperatures studied here. 

The initial kinetic energy of the electrons can be re- 
duced to « 10 mK by using a laser with a bandwidth 
equal to the Fourier transform limit of a 10 ns pulse. One 
may be able to decrease this energy even further by excit- 
ing below the ionization limit. In this case one creates a 
dense gas of highly excited cold Rydberg atoms for which 
many-body interactions can cause a phase transition to 
a plasma-like state |Q. In preliminary experiments we 
have observed the formation of free electrons and ions in 
such a system. This will be the subject of future experi- 
ments. 

The technique to produce ultracold plasmas demon- 
strated in this work is applicable to any atom that can 
be laser-cooled, and other atoms may offer experimental 
advantages. With alkali systems one can attain higher 
initial densities, and alkaline earth ions have accessible 
optical transitions. 

To summarize, we have accessed a new region in the 
parameter space of neutral plasmas by photoionizing a 
cloud of laser cooled atoms. Conditions were realized 
in which both electrons and ions are strongly coupled. 
Experimentally, the initial plasma properties are easily 
controlled and the evolution of the system is described 
with an uncomplicated model. 
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